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We consider magnetization configurations at chiral magnet (CM)/ferromagnet (FM) heterostruc¬ 
tures. In the CM, magnetic skyrmions and spin helices emerge due to the Dzyaloshinskii-Moriya 
interaction, which then penetrate into the adjacent FM. However, because the non-uniform magne¬ 
tization structures are energetically unfavorable in the FM, the penetrated magnetization structures 
are deformed, resulting in exotic three-dimensional configurations, such as skyrmion cones, sideways 
skyrmions, and twisted helices and skyrmions. We discuss the stability of possible magnetization 
configurations at the CM/FM and CM/FM/CM hybrid structures within the framework of the 
variational method, and find that various magnetization configurations appear in the ground state, 
some of which cause nontrivial emergent magnetic field. 
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I. INTRODUCTION 

Magnetic skyrmions are topologically protected non- 
coplanar configurations of magnetization. In contrast 
to vortices and monopoles, they can be embedded in 
a uniform magnetization configuration and behave as 
particle-like objects [1]. Since their first observation in 
the metallic chiral magnet (CM) MnSi by neutron scat¬ 
tering [2] and in (Fe, Co) Si by Lorentz transmission elec¬ 
tron microscopy [3], properties of magnetic skyrmions 
have been extensively investigated [4]. The emergent 
electromagnetic field generated by skyrmionic configu¬ 
rations changes the transport property of the conduction 
electrons, resulting in the topological Hall effect [5-10] 
and the electromagnetic induction [11]. The coupling 
between the magnetization and conduction electrons also 
enables us to control the motion of skyrmions by an elec¬ 
tric current [12-15]: Due to the topological nature of 
skyrmions, they robustly survive in dynamics, and the 
mobility is much higher than that of magnetic domains 
and helical configurations. 

Magnetic skyrmions are observed in CMs, such as 
metallic MnSi [2, 16, 17], Fei_xCoxSi [3, 18], MnGe [19], 
Fei_a,Mn 3 ,Ge [20], and insulating Cu 20 Se 03 [21, 22]. 
These materials have non-centrosymmetric B20-type 
crystal structures, where the relativistic Dzyaloshinskii- 
Moriya (DM) interaction [23, 24] stabilizes a crystalline 
structure of skyrmions, called a skyrmion crystal (SkX), 
under an external magnetic field. Besides the CMs, a lat¬ 
tice of atomic-scale skyrmions is observed in an Fe mono- 
layer on Ir(lll) using the spin-polarized scanning tunnel¬ 
ing microscopy [25], where the four-spin interaction, as 
well as the DM interaction coming from the strong spin- 
orbit coupling in the Ir substrate, plays a crucial role 
to stabilize skyrmions. The enhancement of the DM or 
the spin-orbit interactions at interfaces and, thereby, the 
emergence of atomic-scale skyrmions are actively stud¬ 
ied recently [26-38]. Frustrated spin-exchange interac¬ 


tions [39] and nanopatterned magnetic thin film [40, 41] 
are also predicted to accommodate a stable SkX. 

Mathematically, a magnetic skyrmion is a two- 
dimensional (2D) configuration of a three-dimensional 
(3D) unit vector field (whose manifold is isomorphic to 
the two-sphere 5^), which is classified by the second ho- 
motopy group 7r2(5'^). Hence, the skyrmions observed 
in bulk GMs are cylindrical configurations of skyrmionic 
structures, which are not stabilized in the ground state 
but appear in a small region in the B-T phase diagram 
just below the ferromagnetic phase transition tempera¬ 
ture [2]. It was predicted that the skyrmion crystal state 
can be the ground state of the two-dimensional GMs [42- 
44]. Gorrespondingly, the region of the SkX phase is 
greatly enhanced down to T = 0 in thin films [3, 19]. 
The SkX phase is further stabilized in epitaxial thin films 
due to the magnetic anisotropy [9, 10]. 3D magnetiza¬ 
tion configurations are theoretically investigated in bulk 
and thin films of GMs, where the multi-q configuration 
in the 3D reciprocal space and twisting of the skyrmionic 
structure are discussed [45-51]. The experimental result 
in Ref. [52] suggests that one of the 3D configurations 
may be realized in the bulk MnGe. Appearance of a ma- 
gentic monopole in merging dynamics of two skyrmions 
is also discussed in Ref. [53]. 

In this paper, we theoretically show that by creating a 
hybrid system of a GM and a ferromagnet (FM), various 
3D magnetization configurations appear in the ground 
state. Here, we consider a thin GM and assume that the 
magnetization is uniform along the z direction (the direc¬ 
tion perpendicular to the CM/FM interface) within the 
CM. As in the case of a 2D CM, helical or skyrmionic 
structures appear in the CM at a low magnetic field. 
However, the presence of the FM influences the CM as in¬ 
dicated by the reduced critical magnetic field below which 
skyrmions appear. The non-uniform structures appear¬ 
ing in the CM penetrate into the adjacent FM at a short 
distance from the interface. Hence, a simple helix and a 
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skyrmion-cylinder crystal (SCyX) appear when the FM 
is thin, which are the uniform configurations along the 2 ; 
direction and essentially the same as the spin helix and 
SkX in a 2D CM. As the thickness of the FM increases, 
these structures become unstable and deform inside the 
FM. When a spin helix arises in the CM, the helical 
structure is unwound in the FM by three-dimensionally 
rotating the magnetization vector, forming a sideways 
half-cylinder skyrmion per helical period, which we call 
a sideways-skyrmion array (SSA). On the other hand, for 
the case when a SkX appears in the CM, the skyrmionic 
structure shrinks as one goes deep inside the FM, end¬ 
ing up with a singularity of the magnetization, that is, a 
monopole. We call a crystalline structure of such config¬ 
urations a skyrmion-cone crystal (SCoX). We also con¬ 
sider the case when the FM is sandwiched between two 
CMs with opposite signs of the DM interactions. In this 
case, the helical or skyrmionic structures with opposite 
helicities appearing in the two CMs are continuously con¬ 
nected by twisting the magnetization vector in the FM 
along the z direction, resulting in a twisted helix (TH) 
or a twisted-skyrmion crystal (TSX). 

By minimizing the total energy for each configuration 
mentioned above, we obtain the ground-state phase dia¬ 
grams of the CM/FM and CM/FM/CM hybrid systems. 
Here, we consider within a framework of the variational 
method where the skyrmion radius, the helical pitch, 
and the penetration depth of the non-uniform structure 
are used as variational parameters. We also calculate 
the emergent magnetic field that effectively acts on con¬ 
duction electrons in the strong coupling limit, and find 
that the emergent magnetic field takes nontrivial con¬ 
figurations due to the z dependence of the helical and 
skyrmionic structures: For example, the TH induces a 
staggered emergent magnetic field, whereas the emer¬ 
gent magnetic field for the SCoX points to or from the 
monopole and diverges at the monopole. 

The rest of this paper is organized as follows. In Sec. H, 
we review the phase diagram of a 2D CM with defining 
the characteristic energy and length scales. The vari¬ 
ational method used in the subsequent sections is also 
introduced in this section. In Secs. HI and IV, we dis¬ 
cuss the ground-state magnetization configurations, to¬ 
gether with the emergent magnetic field, at CM/FM 
heterostructures and CM/FM/CM hybrid structures, re¬ 
spectively, by comparing the energies of possible magne¬ 
tization configurations. Discussions and conclusions are 
given in Sec. V. 


II. GROUND-STATE PHASE DIAGRAM OF A 
2D CM 

We first review the ground-state phase diagram of a 
2D CM. We consider a thin film of CM with thickness a 
and assume that the magnetization along the thickness 
direction (the z direction) is uniform. We also assume 
that the emergent structure in the xy plane is much larger 


than the atomic scale so that the energy of the system is 
described using the continuum model as [54] 


E = a j j dxdy 

y (Vn)2 - Dn-{\/ xn) + B{1 - n^) 


(1) 


where n{x, y) is a unit vector describing the direction of 
the local magnetization, J(> 0) and D are the strengths 
of the spin-exchange and DM interactions, respectively, 
and B is the external magnetic field applied in the z 
direction. Here, the origin of the energy is chosen so 
that the ferromagnetic state, n = z, has zero energy. 
The DM interaction favors a non-uniform magnetization 
configurations (V x n| |n), whereas the spin-exchange and 
Zeeman terms are minimized for a uniform magnetization 
aligned in the z direction. The competition between these 
terms results in the nontrivial magnetic structure of SkX. 

A. Spin Helix 

In low magnetic fields, a spin helix appears as a ground 
state, where the magnetization n winds lying in the per¬ 
pendicular plane to the modulation vector k so that it 
satisfies V x njjn. Taking k = kx, the magnetization 
profile is given by 

f 0 \ 

n{x,y) = — sinfcx , (2) 

\ — cos kx j 

whose energy is obtained as 

E^,X2T){k) = aL^ (^k^ - Dk + , (3) 

where L is the system size in the x and y directions. 
Here, the phase of the helix in Eq. (2) is chosen so as to 
satisfy n(0,0) = —z for the sake of convenience in the 
latter discussions. By minimizing Eq (3) with respect to 
fc, we obtain the optimized wave number and energy as 

fchei2D = -y, (4) 

-^hel2D = 

B. Skyrmion Crystal 

Since the helical structure has no net magnetization, 
it cannot survive under a large B and instead, skyrmions 
appear. We first consider an isolated skyrmion. The 
magnetization profile around a charge —1 skyrmion is 
described in the 2D polar coordinates (r, y}) as 

( sin0(r) cos(v3 -I- (j))\ 
sin0(r) sin(v3 -|- />) , (6) 

cos0(r) j 
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(b) 

( 71/^)2 


FIG. 1. (Color online) (a) Magnetization profile of a single 
skyrmion projected onto the xy plane, where the colors on 
the arrows indicate the valnes of riz. The profile is given by 
Eq. ( 6 ) with (j) = — 7 r /2 and 6{r) = 7 r(l — r/^), which has the 
nonzero skyrmion charge, 37 // dxdyn ■ (dxn x dyu) = — 1 . 
The solid circle shows the radius of the skyrmion, r = ^. (b) 
Emergent magnetic field for the magnetization configuration 
shown in (a). 



where 


Aj = 27r / pdp 


Ad = —27r / pdp 


d9o 

dp 


sin^ 9o{p) 


d9o sin 200 (p) 

^ 


2p 


= 27r / pdp[l — cos9q{p)]. 
Jo 


( 8 ) 

(9) 

( 10 ) 


Using 9o{p) = tt{ 1 — p), these coefficients are given by 
Aj = 7 r[ 7 r^ + 7 — Ci(27r) + log(27r)] ^ 7 r( 7 r^ + 2A4:),Ad = 
TT^, and Ab = 7 r(l — d/vr^), with 7 being the Euler’s 
constant and Ci(z) the cosine integral fnnction. From 
the second term in the right-hand side of Eq. (7), we find 
(j) = — 7 r /2 {(j) = 7 r/ 2 ) for U > 0 (U < 0). In the following 
discussions, we choose D > 0 without loss of generality. 

When the energy for a single skyrmion becomes nega¬ 
tive, the system tends to create more skyrmions. Hence, 
skyrmions in the ground state form a crystalline struc¬ 
ture. The total energy for the SkX is evaluated by mul¬ 
tiplying the number of skyrmions to Eq. (7): 


EskxiO 


aL^ /AjJ AdD 


( 11 ) 


where we used (j) = —-k 12. Minimizing Eq. (11) with re¬ 
spect to the optimized skyrmion radius and the energy 
of the SkX are respectively given by 


CskX 


Aj j 

A^D' 


ZT’O _ 

-^SkX — 


aL^ 

TT 


Al,D^ 

2AjJ 


4- AbB 


( 12 ) 

(13) 


where </) is a constant independent of r and (p, and 9{r) is 
a monotonically decreasing function that satisfies 0 ( 0 ) = 
TT and 0 (r) = 0 for r > ^ with f being the radius of the 
skyrmion. The magnetization profile described by Eq. ( 6 ) 
is shown in Fig. 1(a). Here, we introduce a function 
9o{p) (0 < p < 1 ) such that 0 (r) = 0 o(r/^) for 0 < r < ^; 
00 (p) monotonically decreases and satisfies 0 o ( 0 ) = tt and 
0o(l) = 0. In the following discussion, we fix 0o(p) = 
7 r(l — p) and take the skyrmion radius ^ as a variational 
parameter. 

Substituting Eq. ( 6 ) in Eq. (1), the energy for a single 
skyrmion is given by 


For 0o(p) = 7r(l — p), we obtain ^skx — 3.9J/D. On 
the other hand, a SkX is described with a superposition 
of three helical spin textures with the modulation vec¬ 
tors satisfying |fci=i, 2 , 3 | = fchei 2 D and = 0 P], 

which leads to a skyrmion radius (a half of the triangu¬ 
lar lattice constant) ^ = 27r/(i/3fchei2D) — 3.6J/D. The 
small difference between these values suggests that the 
actual profile of 0o(p) does not so deviate from 7r(l — p). 


C. Phase Diagram 

By comparing the energies for the spin helix [Eq. (5)], 
the SkX [Eq. (13)], and the ferromagnetic state [Eferro = 
0 ], the magnetic structure in the ground state changes as 


0 < H < Hcri : 

helix. 

(14) 

Hcrl ^ B Bqj.2 ■ 

SkX, 

(15) 

Bcy2 < B : 

ferromagnet. 

(16) 


EskiiO = “ 



-I- AdD^ sin (j) + AbB^ 



where we have assumed the system size is larger than 
the area of a skyrmion (T^ > tt^s^x); th® critical 




























D. Emergent Magnetic Field 
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ferromagnetic (F) 



skyrmion crystal (SkX) 




One of the striking effects of the appearance of the 
SkX is that it causes the emergent electromagnetic field, 
which then leads to the topological Hall effect and the 
electromagnetic induction [4]. Suppose that the conduc¬ 
tion electron spin is coupled to, and forced to be parallel 
to, the localized magnetization. In the strong coupling 
limit, the electrons behave as if there is an emergent elec¬ 
tromagnetic field defined by 

(-Bem)i = ^ ^ iijkn ■ [djU X dkn) , (19) 

j,k=x,y,z 

{Eem)i = n ■ {din x dtn ), (20) 


where di = djdxi and is the totally antisymmetric 
tensor in three dimensions. In the static magnetization 
configuration, we have £lem = 0, whereas JBem is non¬ 
vanishing for non-coplanar configurations. Indeed, we 
obtain Bem = 0 for the spin helix [Eq. (2)] and 


sin Q dQ ^ 

Bern — 

r dr 


( 7r\ ^ sin(7rr/^) ^ 

?/ W'C 


( 21 ) 


for the skyrmionic configuration [Eq. (6)], where in the 
last equality in Eq. (21) we used 9{r) = 7r(l — r/^). 
The distribution of the emergent magnetic field given by 
Eq. (21) is shown in Fig. 1(b). 


FIG. 2. (Color online) Schematic phase diagram of a two- 
dimensional chiral magnet and magnetization structure in 
each phase. 


magnetic fields are defined as 


Bcrl 

Bcr2 


TT - Ajj/Aj 
tt-Ab 2 J ’ 


AlD^ 

2AbA,jJ 


( > ^crl)- 


(17) 

(18) 


The schematic phase diagram is shown in Fig. 2. Al¬ 
though the actual profile of 0o{p) depends on B, our 
variational function with a fixed 0o{p) can capture the 
ground-state property of the 2D CM. In particular, the 
critical values obtained by using 9o{p) = 7r(l — p) are 
Bcrl = 0.24B^/J and Bcr 2 = 0.67D^/J, which reason¬ 
ably agree with the numerically obtained ones Bcri = 
0.23BVJ and B „2 = 0.780“^/J [14, 55]. 

We should remark here that the above discussion is 
valid only for a thin film as the SkX phase disappears 
from the ground-state phase diagram when a ^ Cskx- 
This is because a conical structure, which is a spin he¬ 
lix along the 2 direction with uniform longitudinal mag¬ 
netization, has lower energy than the SkX. In experi¬ 
ments, the SkX phase is observed in the ground state up 
to a ~ d^skx [22]. 


E. Dimensionless Parameters 


In the following sections, we scale the length in units 
of fskx and the energy in units of AjJL'^/{2TT^skx)- The 
dimensionless variables are denoted with tilde, e.g.. 


X 

Uyi 

II 

(22) 

k = ^Skxfc, 

(23) 

^ 27r^SkxB 

AjJL^ ■ 

(24) 

We also introduce a scaled magnetic field 


-t>cr2 

(25) 

Using these notations, Eqs. (3) and (11) are 
rewritten as 

respectively 


j , (26) 

Bskx(l) = d j + , 

(27) 


and the scaled value for the critical field Bcri is given by 

^ _Bcrl _ AjAb 71" A^/Aj 

Bcr2 TT Ab 


( 28 ) 
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FIG. 3. Schematic of the CM/FM heterostructure, where CM 
and FM denote chiral magnet and ferromagnet, respectively, 
and D is the coupling constant of the Dzyaloshinskii-Moriya 
interaction [See Eq. (29)]. 


III. CM/FM HETEROSTRUCTURE 

We consider a heterostructure of a CM with thickness 
a on a FM with thickness I (Fig. 3). For simplicity, we 
assume that the ferromagnetic interaction is the same in 
the whole system. The total energy of the system is given 

by 


^(Vn)2 + i?(l-n,) 
c? 2 ; yy dxdyn ■ (V x n). (29) 

The 3D magnetization configurations that minimize 
Eq. (29) and the resulting emergent magnetic fields are 
summarized in Figs. 4 and 5, respectively. Here, we as¬ 
sume that the magnetization profile in the CM is uniform 
along the z direction. As we saw in the previous section, 
there are two possible 2D configurations in the CM, i.e., 
the spin helix and the SkX. For each configuration, there 
are two possible configurations in the adjacent FM: When 
the thickness of the FM is small, the magnetization con¬ 
figuration in the CM uniformly penetrates into the FM 
[Fig. 4(a) and 4(b)]; On the other hand, when the FM 
is thick enough, the magnetization configuration is de¬ 
formed in the FM and disappears at a finite depth from 
the CM/FM interface [Fig. 4(c) and 4(d)]. We calculate 
the energy and the emergent magnetic field for the each 
configuration in Secs. IIIA-IIID and discuss the phase 
diagram in Sec. HIE. 




A. Spin Helix 


We first consider the case when a spin helix appearing 
in the CM uniformly penetrates into the FM. When the 
magnetization vector is given by Eq. (2) for all — 1 < z < 
a, the total energy is given by 


Ehel3D(fc) = 


{a + 


aDk -|- (a -f l)B 


, (30) 


or, equivalently. 


Ehei3D(fc) = 7r(a-|-0 



d 2k 
d +1 Ad 



(31) 


By minimizing Eiiei 3 D(fc) with respect to k, we obtain the 
optimized wave number and energy as 


fchel3D 


d Aj 
a 1 Ad 


Q 1 


khel2D, 


^hel3D — Tr{d + l) 


b 

( “ ' 

Aj' 

Ab 

\ a + fy 

Al_ 


(32) 

(33) 


As one can see from the comparison between Eqs. (3) 
and (30), the effective DM interaction relative to the fer¬ 
romagnetic and Zeeman interactions in the CM/FM het¬ 
erostructure is decreased by a factor d/(d -|- /), resulting 
in the reduction of the optimized wave number as shown 
in Eq. (32). As in the case of the 2D CM, there is no 
emergent magnetic field for the spin helix [Fig. 5(a)]. 


B. Skyrmion-Cylinder Crystal 


Similar to Sec. Ill A, when a skyrmion appears in the 
CM and uniformly penetrates into the FM, the magne¬ 
tization vector is given by Eq. (6) for all —l<z<a. A 
possible candidate for the ground state is the crystalline 
structure of such configurations, i.e., the SCyX. The en¬ 
ergy for the SCyX is calculated in the same manner as 
Eq. (11). In the present case, the system size along the 
z direction is elongated by a factor (a -I- ^)/a and the ef¬ 
fective DM interaction relative to the other interactions 
is reduced by a factor a/(a -I- 1). As a result, we obtain 


EsCyxiO 


{a + l)L^ 

TT 


AjJ 

2^2 


a AdD 

(X -\-1 ^ 


AbB 


(34) 


which reduces to 

4cyx(l) = [d + l){l- + b\ . (35) 

VC a + l^ J 


By minimizing Escyx(C) with respect to C, we obtain the 
optimized skyrmion radius and energy as 

IsCyX = ~ ^ , (36) 

/ ~ \ 2 ' 


-^SCyX — {d + l) 


b- 


-l- / 


(37) 


The emergent magnetic field for a single skyrmion 
cylinder is the same as that for a skyrmion in a 2D CM 
and given by Eq. (21). Figure 5(b) shows the configura¬ 
tions of Hem for the SCyX. 
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(c) sideways-skyrmion array (SSA) (d) skyrmion-cone crystal (SCoX) 



FIG. 4. (Color online) Possible three-dimensional magnetization structures in chiral magnet (CM)/ferromagnet (FM) het¬ 
erostructures, where the magnetization vectors in the FM are shown. The magnetization in the CM, which is shown with a 
gray cuboid, is assumed to be uniform along the « direction and forms a spin helix [(a) and (c)] or a skyrmion crystal [(b) and 
(d)]. The magnetization configurations in (a) and (b) are independent of 2 , whereas those in (c) and (d) deform as a function of 
z and become uniform at a distance from the CM/FM interface. In (c), the rows indicated by A and B have opposite whirling 
patterns in the xz plane, whose magnetization profiles are given by n+{x, y, z) and n_ {x, y, z) defined in Eq. (38), respectively. 
The whirling direction is spontaneously chosen for row by row. 


C. Sideways-Skyrmion Array 


When the FM is thick enough, it is not energetically 
favorable to keep the non-uniform magnetization config¬ 
uration in the whole FM. The non-uniform configura¬ 
tion appearing in the CM penetrates only into a finite 
depth of the FM. For the case when a spin helix ap¬ 
pears in the CM, the helical structure is unwound by 
three-dimensionally rotating the magnetization vector as 
shown in Fig. 4(c). This is nothing but a one-dimensional 
array of sideways half-cylinder skyrmions. We pick up 
one of the sideways skyrmion at y = 0 and consider the 


following ansatz: 


n±{x,y,z) 


0 


—sgn(a;) sini?(|a;|) 
cosi?(|a;|) 


(0 < 2 < a). 


= 


'^sin i?( 77 ) cos(±x + 

sin 1 ^( 77 ) sin(±x -|- cj)) | {—I < z <0), 

cos 7 ^( 77 ) 


(38) 


where 77 = + {^z/d)^, x = arg(a; -b i^z/d), (j) = 

— 7 r/ 2 , and 7 ^( 77 ) = 6 * 0 ( 77 /^). Here, we consider an el- 
liptically deformed skyrmion and ^ and d are the ra¬ 
dius in the x and z directions, respectively. We choose 
Oo{p) = 7 r(l — p). Then, Eq. (38) at z > 0 coincides with 
Eq. (2) with k = tt/^. In Eig. 6 , we plot the magneti¬ 
zation vector field given by n^{x, 77 , z) in Eq. (38), from 
which one can see that the helical structure is continu- 
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FIG. 5. (Color online) Emergent magnetic field Bem for the magnetization configurations shown in Fig. 4, where H, SCyX, 
SSA, SCoX stand for helix, skyrmion-cylinder crystal, sideways-skyrmion array, and skyrmion-cone crystal, respectively. Bem 
for (b), (c), and (d) are given by Eqs. (21), (44), and (52), respectively. In each panel, the color scale for the emergent magnetic 
field is scaled by its maximum, where ^ and d are the skyrmion radius (or the half of the helical pitch) and the penetration 
depth of the non-uniform structure, respectively. The direction of the emergent magnetic field is schematically shown with 
white arrows, (a) There is no emergent magnetic field in any direction for the spin helix, (b) For the SCyX, the emergent 
magnetic field in the —2 direction arises at around the center of the skyrmion. (c) For the SSA, the emergent magnetic field 
arises in the y direction and its sign depends on the whirling direction in the xz plane. The rows indicated by A and B in 
(c) corresponds to those in Fig. 4(c), which have opposite whirling patterns in the xz plane, and hence, the direction of the 
emergent magnetic field is opposite, (d) The emergent magnetic field for the SCoX is non-collinear and parallel to e, a unit 
vector along r — with Vm being the position of the nearest monopole. The amplitude of the emergent magnetic field shown 
in (d) diverges at the monopoles. 


ously deformed to a uniform one and that this is indeed 
a half of a skyrmion. Note that though similar configura¬ 
tions have been considered in Refs. [56, 57], the present 
configuration is distinct from them as the axis of the 
skyrmion and, thereby, the emergent magnetic field are 
perpendicular to the external magnetic field in the SSA, 
whereas in Refs. [56, 57] the axis of the skyrmions are 
parallel to the external magnetic field. 


The topological charge for a sideways skyrmion is de¬ 
fined by the integral of the skyrmion density in the FM 


region and calculated as 

1 1 

— / dz / dxn± ■ {dzn± x dxn±) = ±-. (39) 

dTT 2 

As we shall see below, within the framework of the vari¬ 
ational method, the energies for the configurations n+ 
and n_ are degenerate. Hence, which configuration ap¬ 
pears is spontaneously determined for row by row [see 
Fig. 4(c)]. 

By substituting Eq. (38) into Eq. (29), the energy for 


















skyrmions Lj (2^) and its dimensionless value is given by 


Ess Aid, =7r 


aiT^ 2a7r 

Aoi 


ab d f db^ 

~A - 1 -^ ^-^ ^— 7" 

Ab 8 ^ 8d 4 

(41) 


Equation (41) has a minimum with respect to d at 
dssA = , ^ =, 

^1 + 2be 

and the total energy as a function of ^ is given by 
-£^SSA(dsSA,0 


(42) 


2d7r ah V1 + 26^^ 




4^ 


(43) 


Equation (42) shows that d is in the same order as ^ and 
decreases as b increases, which means that the sideways 
skyrmions are compressed to the interface so as to re¬ 
duce the Zeeman energy. In order to compair with other 
comfigurations, we numerically minimize Eq. (43) with 
respect to ^ and obtain the energy of the SSA. 

Since a sideways skyrmion is a skyrmion in the xz 
plane, it induces the emergent magnetic held in the y 
direction. By substituting Eq. (38) in Eq. (19), we ob¬ 
tain 


FIG. 6. (Color online) (a) Magnetization profile of a side¬ 
ways half-cylinder skyrmion given by n+{x,y, z) in Eq. (38). 
Shown are the magnetization vectors projected onto the xz 
plane, (b) The same as (a) but the magnetization configu¬ 
ration at 2 > 0 is replaced so as to form a full skyrmion on 
the xz plane. Namely, n in z > 0 is defined as nxix,y,z) = 
-nAx,y,-z) and nyAx,y,A = ny,zix,y, -z). The solid 
circle shows the region of the skyrmion, (x/^)'^ -|- izjdf' = 1. 
The skyrmion charge for the configuration in (b) is given by 
^ dxdzn ■ {dzTi x d^n) = 1. 


-Bem,± = ri± ■ idzn± x dxn±) 
d^ TTT?/^ 


(44) 


where rj = yj -I- (^z/d)^. The emergent magnetic held 
for the magnetization conhguration in Fig. 4(c) is shown 
in Fig. 5(c), where the rows indicated by A and B in 
Figs. 4(c) and 5(c) correspond to each other. 


D. Skyrmion-Cone Crystal 


a sideways skyrmion is given by 


Essiid,0 =2aL^ 


+ L 


J 
2 

AjJ 


- D 


d_^^_ 
i d 


B 


AbB 


^d 


(40) 


Here, the hrst term in the right-hand side of Eq. (40) 
is the energy of the CM, which is given by Eq. (3) with 
replacing the system size to 2^L, whereas the second 
term comes from the EM. The total energy for the SSA 
is obtained by multiplying the number of half-cylinder 


Similar to Sec. IIIC, when a SkX appears in the CM 
and the adjacent EM is thick enough, the skyrmionic 
structure cannot penetrate into the whole EM [Fig. 4(d)]. 
In Fig. 7, we show the 3D structure developed below 
a skyrmion. We call the structure shown in Fig. 7 
a skyrmion cone. When we see the 2D structure of 
the skyrmion cone perpendicular to the z axis, the 
skyrmionic structure shrinks as one goes deep inside the 
EM and eventually disappears at a finite depth d. Note 
that because a skyrmion is a topologically nontrivial 
structure, it cannot disappear under a continuous defor¬ 
mation. Hence, the skyrmionic configuration ends up 
with a defect of the magnetization, that is, a monopole. 

To give a concrete profile of the magnetization, we con¬ 
sider a skyrmion with radius ^ in the region of 0 < z < a. 
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The total energy for a SCoX is obtained by multiplying 
the number of the skyrmion cones L^/(7r^^) to Eq. (46). 
The dimensionless value is given by 

EsCox{d, i) = jz - ^ d- db + j- -\—Y + Psdb, (49) 

where /3j = Bj/Aj ~ 0.40 and /3s = BbIAb = 1/3. 
Equation (49) has a minimum with respect to d at 

dscox = —r^==^=- (50) 

Similar to Eq. (42), dscoX is in the same order as ^ and 
decreases as b increases, which means that the penetra¬ 
tion depth of skyrmions becomes smaller for larger b. At 
d = dsCoX, the total energy is given by 


-^SCox(l) 


a 

e 


2~a ^ 2Jpj{l + PBh~e) 


(51) 


FIG. 7. (Color online) Magnetization profile in a skyrmion 
cone. The skyrmion radius shrinks as one goes inside the 
FM. At the top of the cone, a monopole emerges. 


whose magnetization vector is given by Eq. (6), and as¬ 
sume that the skyrmion shrinks as a function of z and 
disappears at z = —d < 0. The magnetization profile for 
—d < z < 0 is given by Eq. (6) with replacing 9[r) with 
the following z-dependent function: 


e{r, z) = 00 


^Ui\z\/d) 


(45) 


where /(C) is a monotonically decreasing function satis¬ 
fying /(O) = 1 and /(I) = 0, and ^f{\z\/d) describes the 
skyrmion radius at depth \z\. Substituting Eqs. (6) and 
(45) in Eq. (29), the total energy for a skyrmion cone is 
given by 


13scoi(d,C) =a 


AjJ 


— AdD^ + AbB^" 




(46) 


where the first and second lines of the right-hand side of 
Eq. (46) correspond to the energies for the CM and the 
EM, respectively, and we have defined 


Bj = 27r 




Bb = Ab 



dCfiC). 




(47) 

(48) 


Here, we approximate /(C) = 1 — C and 9o{p) = 7r(l — p). 
Then, the above coefficients are given hy B j = 7r^/2 and 
Bb = 7r/3(l - d/vr^). 


In order to compair with other configurations, we numer¬ 
ically minimize Eq. (51) with respect to C and obtain the 
energy of the SCoX. 

Taking into account the z dependence of 0 and substi¬ 
tuting Eq. (6) in Eq. (19), the emergent magnetic field is 
calculated as 


B 


em 


TT 1 ^ sin[7rr/C(z)] r — Tm 
C(z)_ 7i'r/C(z) z + d ’ 


(52) 


where C(-z) = C/d^l/d) = C(1 + z/d) is the z-dependent 
skyrmion radius, r = (x, y, z), and = (0,0, —d) is the 
position of the monopole. The emergent magnetic field is 
non-vanishing inside the skyrmion cone. It points to the 
monopole, and the amplitude diverges at the monopole. 
The configuration of B^m for the SCoX shown in Fig. 4(d) 
is depicted in Fig. 5(d). 


E. Phase Diagram 

By comparing the energy for each configuration, the 
phase diagram of the CM/FM heterostructure is obtained 
in the (a, b) space as shown in Fig. 8. Here, we calculate 
for (a) I = 0.5, (b) I = 1.1, and (c) I = 1.5. As ex¬ 
pected, for a small I [Fig. 8(a)], only the helix (H) and 
SCyX phases appear. In this case, the phase boundaries 
are analytically obtained by comparing the energies in 
Eqs. (33) and (37) and Ep = 0, and given by 

^F-scyx = { ; (53) 

\a + lj 

^SCyX-H = f ~ ^1) (54) 

Va -I- // 
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FIG. 8. Phase diagram of a CM/FM heterostructure for (a) 
I = 0.5, (b) I = 1.1, and (c) I = 1.5. Here, F, SCyX, H, 
SCoX, and SSA stand for ferromagnetic, skyrmion-cylinder 
crystal, helix, skyrmion-cone crystal, and sideways-skyrmion 
array phases, respectively. The F-SCyX and SCyX-H phase 
boundaries are given by Eqs. (53) and (54), respectively. 
The other phase boundaries are numerically calculated. The 
dashed curve in (c) indicates the F-SSA phase boundary for 
/ > 1 . 


where &i is the critical magnetic field at / = 0 and de¬ 
fined in Eq. (28). Compared with the case for 2D CMs, 
the critical magnetic fields are suppressed by a factor 
[a/(a + 1)]^ due to the reduction of the effective DM in¬ 
teraction. The phase diagram rapidly changes at around 
1 = 1, where the SCoX phase and the SSA phase arise 
between the ferromagnetic (F) and SCyX phases and be¬ 
tween the SCyX and F phases, respectively [Fig. 8(b)]. 
As I increases further, the regions of the SCyX and H 
phases shrink and eventually disappear for / ^ 1. The 
phase boundaries among the F, SCoX, and SSA phases 
do not depend on 1. The dashed curve in Fig. 8(c) shows 
the F-SSA phase boundary for / ^ 1. 



FIG. 9. Schematic of the CM/FM/GM hybrid structure, 
where GM and FM denote chiral magnet and ferromagnet, 
respectively. The signs of the coupling constant D of the 
Dzyaloshinskii-Moriya interaction are opposite for two CMs. 
[See Eq. (55)]. 


IV. CM/FM/CM HYBRID STRUCTURE 


Next, we put another CM on the other side of the FM 
as shown in Fig. 9. We consider the case when the signs 
of the DM interaction in two CMs are opposite. The 
total energy for this hybrid structure is given by 



dz 


dxdy 





— D f dzn • (V X n) -b D 
Jo 



dzn ■ (V X n). 

(55) 


The possible magnetic structures and the resulting 
emergent magnetic field are summarized in Figs. 10 and 
11, respectively. Since the sings of the DM interactions 
are opposite in the two CMs, the helical and skyrmionic 
structures appearing on the top and bottom CMs have 
opposite helicities, which are continuously transformed 
with each other by rotating the transverse magnetiza¬ 
tion by ±7r about the 2 axis. Hence, a TH [Fig. 10(a)] 
and a TSX [Fig. 10(b)] are the possible candidates for 
the configuration in FM with small I, which we discuss 
in Secs. IV A and IV B, respectively. When I becomes 
larger, as in the cases of the CM/FM heterostructure, 
the sideways half-cylinder skyrmions and the skyrmion 
cones come in the FM from the CM/FM interfaces as 
shown in Figs. 10(c) and 10(d), respectively, which we 
discuss in Sec. IV C. The phase diagram is discussed in 
Sec. IV D. 
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(a) twisted helix (TH) 



(b) twisted-skyrmion crystal (TSX) 


(c) sideways-skyrmion array (SSA) 





(d) skyrmion-cone crystal (SCoX) 


FIG. 10. (Color online) Possible three-dimensional magnetization configurations in a CM/FM/CM hybrid structure, where 
CM and FM denote chiral magnet and ferromagnet, respectively. In (a)-(d), the magnetization vectors in the FM are shown. 
The magnetization in the CMs, which are shown with gray cuboids, are assumed to be uniform along the z direction and form 
spin helices [(a) and (c)] or skyrmion crystals [(b) and (d)]. The helicities of the magnetization configurations in the top and 
bottom CMs are opposite. Panel (e) [(f)] shows the magnetization configuration in the bottom CM for (a) and (c) [(b) and 
(d)]. In (a), the helical structure is twisted by —tt about the « axis as one goes from z = 0 to —I and forms a twisted helix. 
Twisting by tt about the 2 axis is also possible. In (b), the skyrmionic structures are twisted by tt (P) or —tt (Q) about the 
2 axis as one goes from 2 = 0 to —1. The direction of the twist is spontaneously chosen for each skyrmion. (c) and (d) are 
the similar structures as Figs. 4 (c) and (d), respectively, where the non-uniform structures (sideways skyrmion and skyrmion 
cone) come in the FM from the both interfaces. In (c), the sideways skyrmions with opposite skyrmion charges (i.e., opposite 
whirling patterns in the xz plane) are degenerate and randomly chosen for row by row. For the case of (c), the rows indicated 
by A (B) have the skyrmion charge 1/2 (—1/2). 
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0 -*/2 ->t 



FIG. 11. (Color online) Emergent magnetic field Bem for the magnetization configurations shown in Figs. 10(a)“10(d), where 
TH, TSX, SSA, SCoX stand for twisted helix, twisted-skyrmion crystal, sideways-skyrmkion array, and skyrmion-cone crystal, 
respectively. Bem for (a), (b), (c), and (d) are given by Eqs. (60), (65), (44), and (52), respectively. In each panel, the color 
coding for the emergent magnetic field is scaled by its maximum, where k, and d are the wave number of the helix, the 
skyrmion radius (or the half of the helical pitch), and the penetration depth of the non-uniform structure, respectively. The 
direction of the emergent magnetic field is schematically shown with white arrows in (a), (c), and (d), whereas that for (b) is 
depicted in the insets, (a) In the TH, because of the twisting along the z direction, the staggered magnetic field arises in the 
y direction, (b) In the TSX, the z component of Bem is the same as that for the SCyX shown in Fig. 5(b), but the x and y 
components, which are absent for the SCyX, arise due to the twisting. The red (blue) color means that the whirling of the x 
and y components is clockwise (anti-clockwise) as depicted in the left (right) inset, whereas the saturation of the color shows 
the amplitude |Bem|. The insets show the conhguration of the emergent magnetic field in the xy plane, where the colors on 
the arrows indicate 6b = arccos[(Bem)z/|Bem|]. Shown are the conhguration for I = 1. The skyrmions indicated by P and Q 
correspond to those in Fig. 10(b). (c) The same as Fig. 5(c). The rows indicated by A and B correspond to those in Fig. 10(c). 
(d) The same as Fig. 5(d). The red (blue) color means that the emergent magnetic held is pointing from (to) the monopole. 
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A. Twisted Helix 


We consider the following magnetic configuration: 
n{x,y,z) 


= 


0 

— sin kx 
, — cos kx I 


{0 < z < a), 


'± sin kx sin(7rz//)' 

— sinfca:cos(7rz//) | (—1 < z <0), 

— cos kx 


(56) 


0 

sin kx 
, — cos kx j 


{—a — I < z < —/). 


In this magnetic prohle, as we go along the 2 direction 
from z = 0 to z = —I, the helical texture is rotated by 
±7r about the z axis so as to continuously connect the 
spin helices with wave vector kx {0 < z < a) and —kx 
(—a — I < z < —/). The total energy for the TH is 
obtained by substituting Eq. (56) in Eq. (55) as 


^TH(fc) = 2Ehel2D(^) +'^1 


1 


b' 

Aj 2Aj 

\i) 

Ab 


= 7r(2a + 1) 



2a 2k 
2a -\-1 


b 

Ab 


^ 2Ajr 

(57) 



FIG. 12. (Color online) Magnetization profile in a twisted 
skyrmion. The magnetization vectors rotate by —tv about the 
2 axis as 2 changes from 0 to — Z. The green thick curves trace 
the positions of the same magnetization directions. 


B. Twisted-Skyrmion Crystal 


Comparing this with Eq. (26), one can immediately see 
that the effective DM interaction relative to the ferro¬ 
magnetic and Zeeman interactions is decreased by a fac¬ 
tor 2aj(2a + /), and the optimized wave number and en¬ 
ergy are respectively given by 


^TH = 


2d 

2d +1 


=fchel2D, 


Eth = 7’‘(2d -I- 1) 


(58) 


' b 

/ 2d ^ 

Aj' 

7r3 

Ab 

\ 25- ~\~ 1 y/ 

Al 

2Ajl 


(59) 


In the case when SkXs appear in the CMs, the 
skyrmions in the top and bottom CMs have opposite he- 
licities, i.e., </) in Eq. (6) is ^ = —tt/ 2 {(j) = (j)/2) for 
0 < 2 < a (—a — I < z < —/). These structures are 
topologically equivalent and can be transformed to each 
other by a continuous transformation. A possible struc¬ 
ture is given by Eq. (6) with taking into account the z 
dependence of (f) as shown in Eig. 12. 

Rewriting (j){z) = g(\z\/l) where g[C,) is a monotoni- 
cally increasing or decreasing function satisfying g(0) = 
—7r/2 or 37r/2 and g(l) = 7r/2, the total energy for the 
crystalline structure of twisted skyrmions is given by 


The last term in the right-hand side of Eq. (59) rep¬ 
resents the additional ferromagnetic interaction energy 
associated with the z dependence of the magnetization 
profile. 

The emergent magnetic field for the TH is calculated 
from Eqs. (56) and (19) as 

Bern = 

where the double sign corresponds to that in Eq. (56). In 
contrast to a simple helix, which has no emergent mag¬ 
netic field, the staggered magnetic field arises for the TH 
as shown in Eig. 11(a) due to the z dependence of the 
magnetization configuration. 


£^tsx(C) =2Eskx(C) + M ^ ^ ^ 


={2h + l) 


1 

F 


2a 2 


2a + I ^ 


11 

I 


(61) 


where 




pdps\ael[p), 


(62) 


and the last term of the most right-hand side of Eq. (61) 
comes from the z derivative of the magnetization. We 
choose 5 (C) = 5+(C) = 7r(C - 1/2) or 5 _(C) = 7r(2/3 - C) 
and 9q{p) = 7r(l — p), obtaining yj = /{2Aj) ^ 0.40. 
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Minimizing Eq. (61) with respect to the optimized ^ 
and the minimum energy are obtained as 




2h -t- I 
2& ’ 


^Tsx — (2a + 0 


b- 


2a 

2a +1 


(63) 


+ y ■ (64) 


As in the case of the previous sections, the DM interac¬ 
tion energy relative to the other interaction energies are 
reduced by a factor 2a/(2a-|-0j ^-nd hence, the skyrmion 
radius becomes larger as I increases. For our choice of 
g{() and 0o{p), the last term in Eq. (64) coincides with 
that in Eq. (59). 

The emergent magnetic field for the TSX is given by 

where the double sign corresponds to that of g±{z). Here, 
the longitudinal component is the same as that of the 
SCyX shown in Fig. 5(b). However, the emergent mag¬ 
netic field of the TSX also has the x and y components, 
which are whirling in the clockwise or anti-clockwise di¬ 
rection depending on the direction of the twisting in the 
EM [see the instes of Fig. 11(b)]. Since the energies for 
the configurations with opposite twisting are degener¬ 
ate, the direction of the twisting is randomly chosen in 
each skyrmion within the framework of the variational 
method, as in the case of the SSA. 


C. Sideways-Skyrmion Array and Skyrmion-Cone 
Crystal 


D. Phase Diagram 


By compareing the energy for each configuration, we 
obtain the phase diagram of the CM/FM/CM hybrid 
system as shown in Fig. 13, where we calculate for (a) 
I < 1.0, (b) I = 2.0, and (c) I = 2.2. When I < 1.0, only 
the TSX, TH, and F phases appear. The F-TH, F-TSX 
and TH-TSX phase boundaries are given by 


6f-th 

^F-TSX 

^TH-TSX 


AbAj 

■^D 


2d 

2d + l 


2d 

2d + l 


2q \ ^ _ IjAb 
2a-\-l) 7il{2d~\~l^ 

2 

1 - 

l{2d +1) ’ 

2 

I &1, 


( 66 ) 

(67) 

( 68 ) 


respectively, where we have used qj = 7r^/(2Mj), and bi 
is defined in Eq. (28). In contrast to the case of Fig. 8(a), 
where the H and SCyX phases start from a = 0, there is a 
lower bound of a for the appearance of the TH and TSX 
phases. For example, From Eq. (66), the F-TH phase 
boundary at 6 = 0 is given by 


( TT / TT^ 2p \ 

" w ^ V ^ ^ 


(69) 


which behaves as 


a 

a 


V2aJ 

t^Ad 

2^f2A.j 


j (F0), 

0.28 {I oo). 


(70) 

(71) 


As in the case of the CM/FM heterostructure, the 
sideways skyrmions and the skyrmion cones may appear 
at the CM/FM interfaces. The resulting structures are 
shown in Figs. 10(c) and 10(d). The energies for these 
structures are twice of those obtained in Secs. HI C and 
HID. 

Since the Bem for a single sideways skyrmion is given 
by Eq. (44), the emergent magnetic field for the SSA 
[Fig. 10(c)] is as shown in Fig. 11(c). For the case of 
the CM/FM/CM hybrid system, the sideways skyrmions 
also appear from the bottom of the FM. The emergent 
magnetic field arises in the +y or —y direction depending 
on the whirling direction of the magnetization vector on 
the xz plane, which is randomly chosen for row by row. 

The emergent magnetic field for the SCoX [Fig. 10(d)] 
is shown in Fig. 11(d). As in the case of Fig. 5(d), the 
emergent magnetic field diverges as one approaches the 
monopole. The crucial difference from Fig. 5(d) is, how¬ 
ever, that the direction of the emergent magnetic field is 
dependent on which interface the skyrmion cone comes 
out from: For the skyrmion cone coming out from the top 
(bottom) CM/FM interface, the emergent magnetic field 
points to (from) the monopole on the top of the cone. 


This is because the ferromagnetic interaction energy as¬ 
sociated with the z dependence of the magnetization 
structure prevents the system from creating non-uniform 
structure. In order to overcome the energy cost in the 
FM, the thickness of the CMs, which have the negative 
DM interaction energy, should be large enough. 

As I increases [Fig. 13(b)], the SCoX and SSA phases 
arise between the F and TSX phases and between the 
TH and TSX phases, respectively, and the regions of the 
TH and TSX phases rapidly shrinks [Fig. 13(c)]. 

V. DISCUSSION AND CONCLUSION 

We have discussed possible magnetization configura¬ 
tions in ground states at CM/FM and CM/FM/CM hy¬ 
brid structures. The energy of the system is calculated 
by using a variational method, where we assume a cer¬ 
tain magnetization structures and take its length scales, 
i.e., the skyrmion radius and the penetration depth, as 
variational parameters. By comparing the obtained en¬ 
ergies, the ground-state phase diagrams of CM/FM and 
CM/FM/CM hybrid structures are obtained as shown in 
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0 12 3 4 

a 

FIG. 13. (Color online) Phase diagram of the magnetization 
structnres shown in Fig. 10 in a CM/FM/CM hybrid struc¬ 
ture with (a) I = 1.0, (b) I = 2.0, and (c) I = 2.2, where 
F, TSX, TH, SCoX, SSA stand for ferromagnetic, twisted- 
skyrmion crystal, twisted helix, skyrmion-cone crystal, and 
sideways-skyrmion array phases, respectively. The F-TH, F- 
TSX and TH-TSX phase bonndaries are given by Eqs. (66), 
(67), and (68), respectively, whereas the other phase bonnd¬ 
aries are numerically calculated. The F-SCoX, F-SSA, and 
SCoX-SSA phase boundaries are the same as those for the 
CM/FM system and independent of 1. The dashed curve in 
(c) indicates the F-SSA phase boundary for / ^ 1. 


Figs. 8 and 13, respectively, where 3D exotic configura¬ 
tions, such as SSA, SCoX, TH, and TSX, appear in low 
magnetic fields. 

In particular, the interface introduces a sort of frus¬ 
tration and hence can produce nontrivial magnetization 
textures absent in each constitute alone. For example, 
both helix and ferromagnet are not topological, while the 
SSA which emerges at the interface between these two is 
topological characterized by the emergent magnetic field. 
The phase diagrams in Figs. 8 and 13 will provide a basis 
to design these nontrivial magnetization structures in the 
interface systems. 

Transport properties of conduction electrons coupled 
to the magnetization is greatly influenced by the emer¬ 
gent electromagnetic field. The distribution of the emer¬ 
gent magnetic field J3em shown in Figs. 5 and 11 will 
produce various topological Hall effects depending on 
the direction of B^m- Especially, the diverging B^m at 
the monopole is expected to affect the electron motion 
strongly. (Note that the lattice constant gives the cut-off 
for this divergence in real systems.) Furthermore, the 
current-driven motion of the magnetization textures via 
the spin transfer torque is determined by the gyro-vector 
G, i.e., the integral of Bem over the space. G enters 
into the Thiele’s equation and the finite G enhances the 
spin transfer torque effect [14]. Once the current-driven 
motion of Bem occurs, the emergent electric field Bem is 
induced, i.e., emergent electromagnetic induction. The 
design of these varieties of phenomena in the heterostruc¬ 
tures will open a rich physics of magnetic textures. 
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